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During cosmological inflation, quasiclassical perturbations are permanently gener- 
ated on super-Hubble spatial scales, their power spectrum being determined by the 
fundamental principles of quantum field theory. By the end of inflation, they serve 
as primeval seeds for structure formation in the universe. At early stages of inflation. 
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■ such perturbations break homogeneity and isotropy of the inflationary background. 



In the present paper, we perturbatively take into account this quasiclassical back- 
ground inhomogeneity of the inflationary universe while considering the evolution of 
small-scale (sub-Hubble) quantum modes. As a result, the power spectrum of primor- 
Q ' dial perturbations develops statistical anisotropy, which can subsequently manifest 

I , itself in the large-scale structure and cosmic microwave background. The statisti- 

, cally anisotropic contribution to the primordial power spectrum is predicted to have 



almost scale-invariant form dominated by a quadrupole. Theoretical expectation of 
the magnitude of this anisotropy depends on the assumptions about the physics in 
the trans-Planckian region of wave numbers. 



. PACS numbers: 98.80. Cq, 98.65.Dx, 98.70.Vc, 98.80.Es 

O 
O 



I. INTRODUCTION 

The scenario of cosmological inflation [1] was designed to address several puzzles of the 
standard cosmology, most of them connected with the property of large-scale homogeneity, 
isotropy and flatness of the observable universe. It was also realized that an inflationary 
epoch can explain the origin of primordial density perturbations [2]. Specifically, quantum 
fluctuations of the inflaton and metric fields, being blown up by inflation, acquire quasi- 
classical features on large (super-Hubble) spatial scales. The predictions of the inflationary 
theory — most notably, the spatial flatness of the universe and adiabatic primordial per- 
turbations with almost scale-invariant spectrum — are surprisingly well verified by modern 
observations of the large-scale structure (LSS) and cosmic microwave background (CMB). 
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Primordial perturbations which serve as primeval seeds for the formation of the observable 
structure are generated at the late stages of inflation (the last 60 or so e- foldings). However, 
depending on a concrete realization of the inflationary scenario, the inflationary stage can be 
rather long, in principle, starting right from the largest classically allowable energy density, 
the Planck density. The total number of inflationary e-foldings can be enormous, with 
typical values ~ 10^^ in the scenario of chaotic inflation. It is clear that generation of super- 
Hubble scalar and tensor perturbations with quasiclassical features cannot be considered 
as specific to the last stage of infiation, but that this process should occur all along the 
infiationary history. This observation is the basis for the theory of the self-regenerating 
eternal inflationary universe [1, 3]. 

Thus one is lead to the idea that the inflationary universe at all times should be in- 
homogeneous on super-Hubble spatial scales. It then makes sense to take this long-wave 
(super-Hubble) quasiclassical inhomogeneity into account when considering the evolution of 
small-scale (sub- Hubble) quantum fluctuations on the background of an inflationary uni- 
verse. Propagating with speed of light, the small-scale quantum fluctuations travel about a 
Hubble distance in a Hubble time, thus passing through inhomogeneities on a Hubble spatial 
scale during inflation. The background inhomogeneity affects the evolution of the small-scale 
modes determining the vacuum state and, as was shown in our previous publications [4], 
the resulting power spectrum of primordial perturbations acquires statistical anisotropy of 
quite speciflc form.^ Namely, the correlation function, say, of the relativistic potential $ at 
the end of inflation is expressed as 



where Pk describes the standard isotropic power spectrum, and i/k is a new anisotropic 
correction, which turns out to be almost scale- invariant. Speciflcally, its multipole expansion 
is given by 



where ns{k) is the spectral index of scalar perturbations (which is known to be close to 
unity), and the quantities A... (A:) depend on k very weakly. Moreover, the expression (2) 
is essentially dominated by the first term, with quadrupole dependence on the unit vector 
n = k/fc. 

In deriving the result (1), only scalar quasiclassical background inhomogeneities during 
infiation were taken into account in [4] . The magnitude of such inhomogeneities in terms of 

^ As noted in [4] , this bears some similarity with the Sachs - Wolfe effect [5] , in which the large-scale inho- 
mogeneity of the universe at a later cosmological epoch is responsible for the locally observed anisotropy 
in the temperature of the cosmic microwave background radiation. 
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the relativistic gravitational potential is estimated to be^ ~ G\H\ (here, G is the Newton's 
gravitational constant, and H is the Hubble parameter). However, inflation also produces 
tensor perturbations [6, 7] with magnitude hi ~ GH"^ on super-Hubble scales, which is 
higher than especially at the early stages of inflation, where \H\ <S H'^. Therefore, 
tensor inhomogeneities of the background of the inflationary universe also necessarily must 
be taken into account, and this will be the subject of the present paper. We will show that 
the results (1) and (2) for the spectrum of primordial scalar perturbations remain to be 
valid; however, in most interesting cases, a major contribution to the quantity i/k comes 
from the tensor inhomogeneities of the background of the inflationary universe. 

We will see that the expected magnitude of the anisotropic part (2) crucially depends 
on the moment of time where one sets the initial (vacuum) conditions for the quantum 
fluctuations of interest at the inflationary state. The earlier this moment is, the larger is 
the number of inhomogeneous quasiclassical modes that have affected the evolution of the 
small-scale quantum fluctuations, hence, the larger is the expected statistical anisotropy. 
This gives rise to an important issue that was first discussed in [4]. It was observed that 
the physical wave number or frequency of a particular quantum mode in the cosmological 
frame exceeds the Planck value Mp ~ during most part of the inflationary universe. A 

natural question arises, whether it is legitimate to use the standard field theory in cosmology 
on spatial scales below the Planck length. Recently, this trans- Planckian issue has become a 
subject of discussion in connection with the possible modification of the form of (statistically 
isotropic) primordial power spectrum [8] and cosmological particle creation [9]. We will 
present several arguments from the literature in favor of the possibility of dealing with the 
trans-Planckian region in cosmology. If the trans-Planckian issue can be safely ignored, 
then, as we will see, the statistical anisotropy can be significant, with the dominant term Aij 
being up to the order of unity. If, however, one is allowed to specify locally homogeneous 
and isotropic initial conditions for quantum modes only after the Planck-radius crossing, 
then the resulting statistical anisotropy is rather small, i/^ ~ 10'^. 

Statistical anisotropy of the power spectrum of primordial scalar perturbations will be 
manifest in the LSS and CMB. The first effect was already discussed in our previous paper 
[4], while the statistical anisotropy of the CMB was a subject of a number of recent studies 
[10-16]. We will briefly describe the arising effects in the present paper. 

Our work represents a theoretical contribution to the issue of the possible origin of statis- 
tical anisotropy in the universe (see [16, 17] for other ideas on this subject). The predictions 
made in this paper have an advantage of being based on the usual theory of chaotic inflation 



^ Throughout this paper, we use the system of units in which h = c= 1. 
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without any additional physical ingredients. Observation of statistical anisotropy in the 
scale-invariant form (2) may thus be regarded as a nontrivial confirmation of the chaotic 
inflation scenario; this is the only possible test of chaotic inflation that we are aware of. 
On the other hand, since the magnitude of the statistical anisotropy depends on both the 
details of the inflationary scenario and the physical assumptions about the trans-Planckian 
region, its nondetection by itself will not invalidate any of the existing models of inflation. 

The paper is organized as follows. In Sec. II, we describe the standard theory of quan- 
tum scalar and tensor perturbations during inflation and introduce our basic notation. In 
Sec. Ill, we discuss the quasiclassical status of perturbations on super-Hubble spatial scales 
during inflation. In Sec. IV, we give a detailed treatment of the evolution of scalar quantum 
fluctuations on small (sub-Hubble) spatial scales on the background of an inflationary uni- 
verse taking into account the effect of scalar and tensor quasiclassical inhomogeneities on 
large (super-Hubble) spatial scales. In Sec. V, we calculate the resulting spectrum of pri- 
mordial scalar perturbations which are seeds for the observable large-scale structure. The 
trans-Planckian issue is discussed in Sec. VI. The observable effects of statistical anisotropy 
of the large-scale structure and cosmic microwave background are described in Sec. VII. The 
results are summarized in Sec VIII. Some calculations are moved to the appendix. 



Throughout this paper, we consider a simple model of inflation based on a single scalar 
fleld </? with canonical Lagrangian and potential V{(fi). As a concrete example which is 
compatible with current observations, we will use the quadratic potential 



with constant mass m^. 

The expanding spatially flat homogeneous and isotropic universe filled with a scalar field 
is described by the set of classical equations 



II. PERTURBATIONS DURING INFLATION 



(3) 



(p + 3H(p + 



dV{cp) 



(4) 



sttg n 



+ vi^) 



(5) 
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H = -47rG02 ^ 



(6) 



where the dot denotes differentiation with respect to the cosmological time t. The so-called 
slow-roll regime of the scalar field is realized when the "friction" term 3H(f dominates over 
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the term with the second derivative in (4), so that we can write 

In view of (5) and (6), the condition for inflation 

\H\ < (8) 

is equivalent to 

0' « V{^) , (9) 
which, for the scalar- field potential (3) during the slow-roll regime (7), reduces to 

A. Scalar perturbations 

In the longitudinal gauge, metric perturbations of scalar type are described by the so- 
called relativistic potentials $ and ^, which enter the metric as follows: 

ds^ = a\T) [{1 + 2^)dT^ - (1 - 2*)dx2] , (11) 

where r is the conformal time. 

Perturbations of the metric are accompanied by those of the inflaton field x = ^V'- For 
perturbations in the linear approximation, one obtains the following system of equations in 
the spatially flat case under consideration: 

* = (12) 
AnGipx ^ + , (13) 



H d 
Hdt 



= 0. (14) 



a? 



aHdt\H 

This system implies a useful relation that expresses $ in terms of x and %': 

(V^ + 4.G^-)$ = l^f^y. (15) 

Here and below, the prime in the time-dependent quantities denotes the derivative with 
respect to the conformal time r, and the symbol is the usual Laplace operator in the 
Euclidean (comoving) three-space. 

Because of the constraint equations (12) and (13), perturbations of scalar type are de- 
scribed by a single scalar degree of freedom. Hence, proceeding to quantization, we note 
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that there exists a single set of creation and annihilation operators realizing the field de- 
composition 



X 



/ 



(27r)3 
(27r)3 



with the commutation relations 



(27r)35(ki-k2) . 



(16) 
(17) 

(18) 



The spatial dependence of the mode functions is given by the standard exponent e*"*^, 
and their temporal behavior, described by the functions Xki^) ^ki^), depends only on 
the absolute value of k, which implies the spatial homogeneity and isotropy of the vacuum 
state defined by the annihilation operators. 

To determine the temporal dependence of the mode functions, hence, the vacuum state, 
one can introduce the variable 

v^a(^X + (19) 

that brings the quadratic action for scalar perturbations into a canonical form^ (after the 
constraints are taken into account) 



.'2 



divd^v H v' 

z 



drd^x , 



(20) 



where z — a(p/H. In the high-frequency regime, where k ^ aH, one deals with the quantum 
theory of a free massless scalar field in Minkowski space with a well-defined quantum state 
determined by positive-frequency solutions. 

For our purposes, however, it will be more convenient to use the variable x the high- 
frequency region. Using constraint (15), one can see from (19) that ~ ax in this region, so 
that the two variables are simply related. One can describe this fact as insignificance of the 
self-gravity of the infiaton field in the high-frequency region. In this case, the appropriately 
normalized positive-frequency modes for the infiaton field have the form 

1 



Xk 



(21) 



aV2k 

The expression for the modes of the relativistic potential in the high-frequency region 
are then obtained using equation (15): 



inGcp ( - 



3/2 



k:$> aH . 



(22) 



Here and below, the spatial indices are raised and lowered with the use of the Euclidean metric 6ij . 
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Here we used the fact that the value of changes slowly with conformal time relative to the 
frequency k ^ aH (during the slow-roll regime, it varies slowly with time even relative to 
H). 

A more careful analysis of equation (14) shows that the approximate solution (22) is, 
in fact, valid in a somewhat broader range of wave numbers, namely, for k ^ a\H\^^'^ (we 
remember that ^ during inflation). Indeed, after introducing a new variable 

0=|ij|-V2$^ (23) 
we can write equation (14) in terms of the conformal time as 

0"+ (F- V')0 = O, (24) 

in which the time-dependent term 



a' 



H HH 3 ( 



2H 2H 4:\H 



+ 2H 



(25) 



is of the order of a^\H\ during slow-roll inflation. Taking into account equation (6), one can 
see that (22) is a good approximate solution of (23), (24) in the domain k » a|ijp/^. Then, 
using (13) and the slow- roll condition, we can write the approximate solutions for the modes 
valid in a broader frequency range as follows: 

/qX 3/2 

f« IttGc^ f - j e-'*^% k:^a\H\^/^, (26) 

X. « ^ f 1 + ^) e-''^ , k » a\H\'/' . (27) 
aV2k V k J 

The Fourier amplitude of the solution of equation (14) can also conveniently be written 

as 

= Ak- / —^dt + Bk- , 28 

where and are some time-dependent functions. As the physical wavelength X — a/k 
of the perturbation increases, the functions and Bk approach constant values. This is 
clear from equation (14), in which one can neglect the Laplacian for k <^ aH. In this limit, 
solution (28) with and Bk which are constant in time becomes exact. 

To write an approximate solution in the region k <S aH, we note that the second term 
in (28) rapidly decreases during inflation and becomes insignificant. The whole expression 
(28) can then be approximated as 

~ Ak— , k^aH. (29) 
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Then equations (12) and (13) imply 

Xfe~-A|:, k<^aH. (30) 

One can determine the amphtude by noting that the frequency domains in (26), (27), 
(29), and (30) overlap. One obtains 

' (31) 



V2A;3/2 ip ' ' ' a 

During the future evolution, the value of should tend to a constant, which is then 
estimated by the value of the ratio H'^/(p at the moment of time where^ a|ij|^/^ = k. We 
will supply the homogeneous fields evaluated at this moment of time by the index "A;" so 
that, for instance, (fk — 0la|ij|i/2=fc- Thus, we have 

P, = k^\Akf^^, a\H\^'''>k. (32) 

In view of definition (29), the last expression gives the primordial power spectrum of 
perturbations for the variable $ at the end of inflation, which is conventionally defined as 
the moment of time where = H^- The quantity typically is a logarithmic function of 
the wave number k. In the slow-roll approximation with potential (3), we have^ 

12^"'^''"^*' H^^^mlip\ (33) 

so that 

a^af e-2-G(^^-^?) ^ (34) 

where the index "/" refers to the moment of the end of inflation, where \H\ — H'^. Combining 
these expressions, we eventually get 

n« f Gmj(i + In ^)'. (35) 

in which 

In ^ 50 (36) 



^ Since the universe expands exponentially rapidly during inflation, the value of the slowly evolving quantity 
iJ^/tp does not change much in the time interval during which < k/a < H. For this reason, the 

amplitude of perturbations is frequently estimated at the moment of the Hubble-radius crossing during 
inflation, i.e., at aH = k. 

^ To our knowledge, relations (33) and (34) were first obtained in [18], where the dynamics of a universe 
dominated by a massive scalar field was under consideration. 
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for the typical wavelengths corresponding to the large-scale structure. Since Phjl'K^ ~ 2.4 x 
10~^ according to the observations of the angular power spectrum of the CMB temperature 
anisotropy [19], we get the following estimate for the mass of the inflaton field: 

Gml ^2x 10"^^ . (37) 

The running spectral index ns of scalar perturbations is conventionally defined as 

so that a scale-invariant power spectrum corresponds to ns = 1. For spectrum (35) and for 
(36), we calculate 

715 1 - 2 ^In ~ 0-96 ■ (39) 

B. Tensor perturbations 

The idea of parametric amplification of quantum tensor fluctuations (gravitational waves) 
in a homogeneous and isotropic universe was proposed in [20] . The spectrum of gravitational 
waves produced during the de Sitter stage in the early universe was first calculated in [6] and 
later analyzed in [7] in the context of the theory of cosmological infiation. Here we provide 
the relevant relations that will be used in the following sections. 

Tensor perturbations can be presented in the form 

where r]^i, is the Minkowski metric, and h^^ denotes the small perturbation. One can work 
in the transverse traceless gauge 

/ioo = h^i = 5'=Kj - , d%j - . (41) 

For a plane wave with wave vector k, there exist two independent polarizations, which are 
usually denoted by cr = -|-, x. 

The classical dynamics of tensor perturbations is determined by the expansion of the 
Hilbert-Einstein action to quadratic order in hij : 

Sg = ^ j a'{r)rj'^''d,0d,Kj drd'^ . (42) 
One can proceed to the Fourier representation 

h^A^,T)= / tS^ E V.(^)e^.(k,a)e^'^^ (43) 
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where eij (k, a) are symmetric traceless polarization tensors with real components satisfying 
the conditions 

Cij (k, cr) = eij{-k, a) , k'-Cij (k, cr) = , e'^ (k, a) (k, a') = 6^^' ■ (44) 
The action then takes the form 

Quantization is performed by interpreting hk^^iT) as operators with the decomposition 

/ik,a(r) = hk{T)a^,, + /i^(T)alk,, (46) 

and commutation relations 

«k,.,, = (27r)^5.,.,5 (ki - k^) . (47) 

A new, canonically normalized, variable 



satisfies the equation 

v'^ + u;1{t)v, ^ , (49) 

where 

^ fc2 _ ^ ^ ^2 _ ^2 ^ 2i/2^ . (50) 

The initial conditions for the normalized solution with positive frequency at the moment of 
time Tj are 

Note that these initial conditions make sense only for sufficiently small wavelengths, for 
which ul > 0, or k'^ > c? [ii + 2if^^ , at the initial moment of time. 

In the inflationary regime, where ^ . the normalized solution which has positive 
frequency in the asymptotic past is given by the approximate expression [compare with (27)] 



v^^^U^'^^ e--^^^ , k » a|ij| V2 . (52) 
V2k V k ) 

It is exact in the de Sitter space, where H = const. On small spatial scales, the nondecaying 
solution of (49), (50) is Vk oc a, so that 
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III. QUANTUM AND CLASSICAL PERTURBATIONS 

Observing expressions (26), (27) and (52), one can notice a qualitative difference between 
the behavior of the vacuum modes on sub-Hubble spatial scales [k ^ aH) and their behavior 
on super- Hubble spatial scales {k <^ aH). On sub-Hubble spatial scales, they are quite 
similar to the corresponding vacuum modes in the Minkowski space: they oscillate with 
positive frequency in time, and the scalar-field amplitude has the characteristic behavior 
oc The field's quantum state in this domain of wave numbers features the usual 

vacuum in a fiat space. On super-Hubble spatial scales, the temporal evolution is "frozen" 
and the scalar- field amplitude contains an extra factor k~^. The property of the quantum 
state allows one to speak about the generation of a quasiclassical field condensate on large 
spatial scales (see [1]). Indeed, the equal-time correlation function, say, for the perturbations 
of the infiaton field, is given by 

(X(xjx(yj) J ^2^^3 IX.I J (2^)3 j,^^ + 2klJ ' ^'^^ 

where 

kph = - , Xph = ax (55) 

are the physical wave number and physical distance, respectively. The first term in the 
brackets of (54) corresponds to the usual quantum fiuctuations in the fiat space. The second 
term, from the viewpoint of fiat space, can be interpreted as the presence of particles with 
occupation numbers 



^ ' ^^^^ 

leading to nonzero correlation in (54) for exponentially large physical distances H <^ 
I Xph — Yphl ^ Ha/ai, where Oj is the scale factor at the beginning of infiation. 

The behavior of the quasiclassical part of the infiaton (/? in the infiationary universe can 
be described in terms of Brownian motion induced by conversion of quantum fiuctuations of 
the field </? into quasiclassical inhomogeneities. The probability distribution of the local value 
of the infiaton field obeys an effective diffusion equation (see [1, 3, 4, 21] for the development 
of this theory) . 

At the end of infiation, perturbations on super-Hubble spatial scales constitute the pri- 
mordial density perturbations. It is usually assumed that they represent a realization of a 
classical Gaussian random field described by the correlation function (54) and accompanied 
by the corresponding classical perturbations in the relativistic potential The mechanism 
of this transition — from the quantum state that does not break the symmetries of the 
background space to a classical state that obviously breaks these symmetries — remains 
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to be an open issue. It is presumably related to the problem of collapse of wave function, 
hence, to the deep interpretational issues of quantum mechanics (see [22]). 

IV. SCALAR PERTURBATIONS ON THE INHOMOGENEOUS 
INFLATIONARY BACKGROUND 

A. Self-regenerating inflation 

It is clear that the generation of scalar and tensor perturbations on super-Hubble spa- 
tial scales that subsequently behave quasiclassically is not restricted specifically to the last 
stage of inflation, where it is supposed to produce the primordial seeds for the subsequent 
formation of large-scale structure and primordial gravity waves, but that this process occurs 
all along the inflationary history. Therefore, the quasiclassical parts of $ and hij are per- 
manently generated during inflation constituting a (random) inhomogeneity and anisotropy 
of the inflationary universe. It thus makes sense to make the next step in the perturbation 
theory and consider the evolution of short-wave {k ^ aH) perturbations (to become the 
observable large-scale structure) on the background of the long-wave {k < aH) quasiclas- 
sical inhomogeneities permanently generated during inflation. In order to set up a valid 
theoretical background for this procedure, we need to address the issue of the evolution of 
the inhomogeneous inflationary universe. This is the subject of the present subsection. 

Because of the permanent generation of classical long- wave perturbations and their subse- 
quent evolution, the inflationary universe considered as a whole may be very inhomogeneous, 
its different spatially remote parts having quite different histories. Indeed, the dispersion of 
the inhomogeneity of the classical inflaton field generated during one Hubble time is equal 
to (x^) ~ if^/27r according to (54). At early stages of inflation, a typical random change 
of the inflaton can be larger than its regular change A(p — (pH~^ caused by slow rolling; a 
criterion for this property is A(p < ^^(x^ , or 



dip 

For the scalar- field potential (3), this condition becomes 



which is considerably higher than the inflationary boundary (10) but still lower than the 
boundary of the Planckian energy density pp ~ in view of (37): 

y((^,) ~ m^G-^/^ « G-^ . (59) 
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If the inflaton field is above the threshold given by (57) or (58), then its value has a 
significant chance to increase with time due to the mechanism of quantum fluctuations 
described above. This random behavior of the classical part of the inflaton field leads to a 
picture of eternal self-regenerating inflation [1, 3]. One of the consequences of this picture is 
that the inflationary universe as a whole cannot be split into a homogeneous and isotropic 
background and a small classical perturbation because very soon one has to deal with 
perturbations with relative amplitudes larger than unity, which are impossible to treat in a 
perturbative way. 

To overcome this difficulty, we will follow the history of a spatially bounded domain of 
the inflationary universe choosing its variable comoving size to be conveniently larger than 
the size of the horizon at every moment of time and, on the other hand, not very large, so 
that inhomogeneity within this domain can still be treated as small. This approach was also 
adopted in [4] . Specifically, in terms of the original coordinates (r, x) , at every moment of 
time, we deal with a domain of comoving size Ax — {(aH)~^, where the constant <^ 1, 
i.e., a domain well exceeding the Hubble size. Within this domain, at every moment of time, 
we can split the quasiclassical perturbations generated during inflation into the spatially 
homogeneous parts $ and hij and the remaining parts with zero spatial average, which we 
denote by $c a-nd /i? (endowing them with the subscript or superscript "c"). After that, 
we eliminate the homogeneous parts l» and hij, depending only on time, by making the 
transformation of coordinates and of the scale factor within the domain under consideration 
as follows. 

To ehminate l», it suffices to proceed to a new time coordinate: — > (1 — 2l»)dr. This 
will redefine the scale factor as a — > a(l — $). After this transformation, the homogeneous 
part of the metric takes the form 

ds^ = [dr^ — (5jj — /ijj) dx^dx^] , (60) 

where hij depend only on r. 

To get rid of the global anisotropy in the spatial part of the metric (60), we make the 
coordinate transformation 

X' X' + \^KjX^ , T ^ T - ^^^jX-'x^ , (61) 

where we set the origin of coordinates x^ to be at the center of the domain under consider- 
ation. This transformation brings the metric (60) into the form 

ds^ = [{1 + 2$;,) dr^ - (1 - 2^;,) cix^] , (62) 

where 

= -| m + ahlj) x'x^ , = ^ah[jx'x^ , (63) 
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and a = a' /a. 

The complete classical metric within the domain under consideration is then given by 



ds^ = [(1 + 2$ft + 2$c) dr^ - (1 - 2^,, - 2^c)dx^ + hljdx'dx^] 



(64) 



Since the quantities $c and describe inhomogeneities with zero spatial average over the 
specified domain, they can be approximated by the following expressions with conventionally 
chosen time-dependent exponential momentum cutoff 



ikx _ ^-k/CaH^^ 



-k/saH 



(65) 



X 



ikx 



•-k.o- 



-k/eaH 



(66) 



where the exponent with the constant £ 1 determines the momentum integration region in 
(65) and (66) ensuring that we are well in the quasiclassical domain, and the exponent with 
the constant C ^ 1 ensures the property of negligible spatial average over our monitored 
spatial domain with comoving size {C,aH)~^. The quantity hij entering (63) is the average 
of the classical part of the corresponding field over the spatial domain under consideration 
and, therefore, is given by the expression 



hi 



{2ny 



J2 %(k,(7) hf^ir)a^,^ + hl{r)al^ 



-k/CaH-k/saH 



(67) 



<T=+,X 



in which the effective momentum cutoff parameter in the exponent is — ( ^ + e ^ 
in view of the condition ( <^ e <1. 



B. Small-scale perturbations on the inhomogeneous inflationary background 

Now we are going to investigate the propagation of small-scale scalar modes on the back- 
ground of the metric (64). In our treatment [4], this was done by taking into consideration 
only scalar classical inhomogeneities $c present in (64). In this paper, we are going to take 
into account both scalar and tensor classical background inhomogeneities. 

For an inhomogeneous background, the convenient classification of perturbations into 
scalar, vector, and tensor types is no longer exact. In particular, if one wishes to use it, 
then one has to deal with mixing and interaction between these types of modes mediated by 
the background inhomogeneity. However, for the modes with extremely short wavelengths, 
much below the Hubble scale, this mixing can be neglected. Indeed, the self-gravity of the 
scalar infiaton field is negligible at such short wavelengths, and perturbations in a scalar 
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field are described only by scalar quantities. We, therefore, use the same method as in [4]. 
Specifically, for the scalar modes in the high-frequency regime (with physical wavelengths 
A <^ H~^), we neglect the effects of self-gravity and consider propagation of the inflaton 
perturbations on the background of the metric g"^^ defined in (64). The corresponding 
small-scale perturbations x = ^f a-re then described by the equation 

1 



where 



--d, (^/^^r^^x) + mjx = , (68) 



mli^^) - ^ (69) 



is the scalar-field effective mass, and (pc denotes the quasiclassical part of the scalar field, 
defined similarly to (65). For the free scalar field, m^{ipi^ will be just the constant inflaton 
mass m^. In what follows, however, we will sec that the mass term can be neglected 
altogether for the high-frequency modes described by (68). 

We have the right to treat equation (68) only in linear approximation in the quasiclassical 
inhomogeneities <l>c and h^^ since our quasiclassical metric (64) is valid only in linear 

approximation in these flelds. Equation (68), written explicitly to linear order in terms of 
the background classical inhomogeneities, reads 

X" - V\ + 2ax' + a'rnjx = (4$;, + + 3^',) x' + (4$^ + 2$, + 2^,) V\ 

+ ($, - *,) d,x + h^^d.d^x] - 2a2mJ$,x , (70) 

where a = a! ja^ and we remind the reader that the spatial indices are raised and lowered 
using the Euclidean metric bij. The mass function (69) for an arbitrary inflaton potential 
is also supposed to be expanded in (70) in the linear perturbation of the inflaton. However, 
as we already mentioned, the corresponding terms in (70) will be insignificant. 
The quantum field has the decomposition quite similar to (16): 



X 



(27r) 



3 



ikx I * — ikx t 

XkC Ok + XkS Ok 



(71) 



with the creation and annihilation operators satisfying the commutation relations (18), and 
the modes Xk(T, x), which now depend on the whole vector k as well as on x, satisfying the 
usual normalization conditions 

(Xk„ XkJ = / i (xL.Xk, - Xk.Xk'J e^^^^-'^^^^y^^^r^'x = (27r)^5 (k^ - k^) . (72) 

J r=const 

Following [4], we look for the solution of (70) in the form 

Xk = Xfce-'^''~''^N (73) 
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where Xki^) solutions (21) for the unperturbed modes, and the (complex) phases Sk{T, x) 
and 7k (t, x) appear due to the influence of the scalar and tensor quasiclassical inhomo- 
geneities, respectively. The dependence of on and of Tk on /i? , will be considered in 
linear approximation. Furthermore, since we are working in the domain of high frequencies, 
we will use the following formal expansion in powers of : 

Su = kSo + iSi + {k-^) , Tk - kTo + iTi + (A;"^) , (74) 

where we have omitted the implicit index k in the functions 5*0, 5*1, Tq, and Ti. In fact, 
this is the expansion in powers of the dimensionless parameter {ki)~^ <^ 1, where £ is the 
comoving size of the quasiclassical background inhomogeneities. With this notation, all the 
functions ^o, Si, Tq, and Ti turn to be real, as we will see soon. 

Substituting (73) and (74) into (70), we get the following system of equations for the 
phases: 

{dr + n- V) 5-0 = 2$e , (75) 
{dr + n- V) Si = 2$', - 1 {S'^ - V'So) , (76) 
{dr + n-V)To = (^^H + '^h + Q, (77) 
{dr + n.V)Ti = 1 + 3*;) - in • V (*,-*,)- ^ (T^' - V^To) , (78) 

where n = k/Zc is the unit vector with components rii, and 

Q = ^hi^n,nj . (79) 
The solution for the components of the phases -S'k [$c] and Tk [h^j] is straightforward: 
5o(t,x) = f 2$e[T,x-n(T-f)]cif, (80) 
Si{t, x) = $e(r, x) + [V^ - (n • V)'] [ $e [f , x - n (r - f)] (r - f) df , (81) 
To{t, x) = /""($, + + g) [f , X - n (r - f)] df , (82) 
Ti{t, x) = n • V g [f, X - n (t - f)] - ^Q{t, x) + ^^(t, x) 

+ i[v2-(n-V)'] y""($;, + *^ + g)[f,x-n(T-f)](T-f)df. 

(83) 



The moment of time Ti in (80)-(83) is taken to be the time at which the initial conditions 
for perturbations are set. It may or may not coincide with the moment of the beginning of 
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inflation (we wiU discuss this issue in Sec. VI). For simplicity, we can assume the background 
quasiclassical inhomogeneities to be absent at this moment of time, i.e., we can set ^c{Ti, x) = 
and h^j{Ti,x) = 0. This enables one to speak about the homogeneous and isotropic 
initial conditions for quantum fluctuations in a homogeneous and isotropic initial universe, 
hence, to ensure that the inhomogeneity and anisotropy are generated intrinsically in a self- 
consistent way. Alternatively, one can consider a region of the initial universe with size much 
smaller than the local Hubble radius and much larger than the wavelength of the small-scale 
modes under consideration, and by making the local change of time and space coordinates 
in that region in the neighborhood of the initial moment of time, one can make sure that 
solutions (80)-(83) with the initial conditions 5'k(Ti,x) = Tk(Tj,x) = correspond to the 
locally homogeneous and isotropic vacuum. (See [4] for more details in the case of scalar 
background inhomogeneity.) 

We trace the evolution of the mode Xk the form (73) until some later moment of time 
which is close (say, a dozen e-foldings) to the moment of its Hubble-radius crossing. At this 
time, the mode is still in the high-frequency regime, but the classical inhomogeneities $c 
and h^j that substantially influenced its evolution during the preceding inflationary period 
are already stretched well beyond the Hubble radius because of inflation. Therefore, one can 
safely neglect the influence of $c and /i? at this stage and return to the purely hnear theory. 
We also need to take into account the self-gravity of the inflaton perturbations since it is 
going to become important after the Hubble-radius crossing. Thus, we need to determine 
the behavior of the mode $k at this period of time, well before the Hubble-radius crossing. 

The modes $k enter the held decomposition similarly to (17): 

_fk_ 
(2^ 

We look for the solution of $k in the form similar to (73): 

= ^^e-'^''-'^'' , (85) 

where ^^(t) is the unperturbed solution (22), and 5'k(T, x) and 7k(T, x) are the phases corre- 
sponding to 5'k and Tk in (73), describing the influence of the scalar and tensor quasiclassical 
inhomogeneities, respectively. We can expand them in powers of similarly to (74): 

Sy, = kSo + iSi + (k-^) , fk = kfo + ifi + {k-^) . (86) 

The equation of the linear perturbation theory that connects <l>k and Xk is (15). We 
substitute (73), (74) and (85), (86) into this equation and, by collecting equal powers of k~^, 
obtain the following expressions for the new phases: 

So^So, Si^Si + S'o + 2n-VSo^Si + n- WSo + 2$^ , (87) 
To = To , fi = Ti + + 2n • VTo = 7^1 + n • VTo + + ^'ft + g . (88) 



*ke^'X + *ke-''M ■ (84) 
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On the next step, we should hnk the solution (85) for the mode well before the Hubble- 
radius crossing to the solution well after the Hubble-radius crossing. The solution in this 
last regime can be expressed similarly to (29): 

*k~^k(x)-^, (89) 

with the quantity >lk(x) depending on x as well as on k. Comparing this equation to (85), 
one can find a nontrivial phase factor in : 

^k(x) = AkC-'^^^^^'"'^-'^'^^^^'''^ , (90) 

which we can evaluate at the moment r^, defined by the condition a\H\^^^ — k, because, as 
we have seen in Sec. II A, the high-frequency approximation for $k can be extended to the 
domain k > a\H\^^'^ [see (26)].^ The isotropic quantity was also defined in Sec. IIA. 
Therefore, at the end of inflation, which we defined by the condition \H\ — if^, we have 

$k (t/, x) ^ -Akc-'^^^^^'"'^-'^'^^^^'''^ . (91) 

This result enables us to calculate the power spectrum of the primordial perturbations. 



V. POWER SPECTRUM OF THE PRIMORDIAL PERTURBATIONS 



A. General expression 



The correlation function of the primordial scalar perturbations ($(x)$(y)) is associ- 
ated with the symmetrized quantum-mechanical average of the corresponding field operator. 
However, in the theory under consideration, the quantum operators Oj^ and a]^ enter the field 
$ not only explicitly in the decomposition (84) but also imphcitly in the phases 3]^ and Tk in 
(85). These phases depend on the field inhomogeneities with wavelengths greatly exceeding 
the size of the observable universe; therefore, they cannot be included in the quantum av- 
eraging procedure. The uncertainty in the correlation function due to the presence of these 
modes is to be regarded as an instance of cosmic variance. 

Using the new expression (91) for the modes, we can write the correlation function at the 
end of inflation 

(*(x)*(y)) = I J -^slAkl'e^^-^^ 

X L-i^k{-rfc.x)-i^l(-^fc.'')+i^i;(Tfc,y)+ifj;(Tfc,y) _^ g-i5_k(rfc,y)-if_k(rfc,y)+i5* k(rfc,x)+if *k(T;t,x)j /ggN 



^ The dependence of the phases in (90) on in the neighborhood of the Hubble-radius crossing is very 
weak, as will become clear later, so the precise choice of the moment of time Tk is not very important. 
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which we symmetrized with respect to the x <-> y interchange. 

Since (Sk and Tk depend on the modes with wavelengths much exceeding the size of the 
observable universe, we can expand these quantities in power series in the coordinate distance 
X in the neighborhood of the observer. Expansions (86), as we said, are, in fact, expansions 
in powers of {M)~^, where I is the comoving scale of inhomogeneity of the phases <5k and 31; 
therefore, we can restrict ourselves to terms at most linear in x in the expansion of S'o and 
To and retain only the zero-order terms in the expansion of Si and T\ (i.e., evaluate them 
at X = 0, which we take to be the center of the region of observation). This procedure will 
produce terms hnear in the difference x — y in the exponents of (92), so that the coordinate- 
dependent parts of the phases will take the form (k — (5k) (x — y). It is then convenient 
to change the integration variables from k to k — (5k in (92). After this cumbersome but 
straightforward procedure (see Appendix A for details) , we arrive at the following expression 
for the correlation function, which we now write in linear approximation in $c ^-nd /i? : 

($(x)$(y)) = j (1 + -k" + -k^) e^^(--) . (93) 

The quantities and z/^ ^'^^ ^^e resulting effects of 5'k and Tk in (92), i.e., of the scalar 
and tensor quasiclassical background inhomogeneities, respectively. They are given by the 
expressions (see Appendix A): 

-\{ns-l)^-V [(5o)k - (^o)-k] , (94) 
= ^ (^5 - 1) n • V [(ro)k - (To)-k] . (95) 

Here, we have restored the index k in the components 5*0 = Sq and Tq = Tq and used the 
definition (38) of the running spectral index. Expressions (94) and (95) are to be evaluated 
at r = Tjfc and x = 0. 



B. The form and magnitude of statistical anisotropy 

1. The scalar part, 

The contribution of the scalar inhomogeneity to the statistical anisotropy, z/^, was de- 
scribed and studied in the previous papers [4]. Here, this analysis will be carried out in 
more detail. Using (65), (80) and (94), we have 

vi^{ns-l)I^, (96) 
= 2 r dr I i^$^(r)(pn) sin [pn (r^ - r)] e'^l^'^^a^ + H.c. , (97) 
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where n — k/k, the scalar product of p and n is denoted by pn, and "H.c." denotes the 
Hermitian conjugate of the preceding expression. 

The variance of is estimated in Appendix B, with the result 

Hence, 

The numerical value of this estimate is obtained by using (39) and setting £ ~ 1: 

'{ly^Y^ ^5x IQ-^G {Hf - Hi) ~ 5 x IQ-^GHf , (100) 

where the last equality assumes Hi » Hk- Thus, if we set our initial conditions at the 
largest possible energy density — the Planck energy density (assuming that inflation starts 
from this energy density), — we will have GHf ~ a few and ~ 10~^. The results (99) 
and (100) were first obtained in [4]. 

Concerning the form of the anisotropic part z/f , one can sec from (96), (97) that it strongly 
depends on the direction n = k/A; but very weakly on the absolute value of k for given n. 
Indeed, the spectral index ng depends on /c in a logarithmic way, as can be seen from (39). 
The derivative of integral (97) with respect to In A; for fixed n = k/A; is given by 



dhik 



(101) 



The variance of this quantity can be estimated quite similarly to the procedure (B3) of 
Appendix B using the relation 

dn _ \Hk\'/' 



dink kHk 

(which one can easily verify), and we obtain the following result: 



(102) 



oij^^ \ \ e 



dink J I IStt 
The ratio of the variances is given by 



G\Hk\ . (103) 



Qjs N 2^ 



d\nk 



< 1 . (104) 



This establishes very weak dependence of on the absolute value of k for fixed n = k/A;. 
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Let us now consider the dependence of the anisotropic part on the direction n — ii/k. 
To this end, we can expand the sine under the integral in (97) in powers of its argument, 
thus obtaining the formal series 

oo 
m=0 

where 

liiK n) = ^1=^ £ dr in - rf-^' J |^$p(r)(pn)^(-+^)e-/-^ap + H.c. (106) 

This is nothing but a multipole expansion in powers of n. 

Proceeding exactly like in the derivation of (B3), we can estimate the variance of every 
partial term in (105): 



2^ 

{&) ^-^^ 2 • (107) 

^ / (4m + 5)(m + l)2[(2m + l)f ^ ^ 

One can sec that the quantity is dominated by the zeroth term Iq , since the next term 
already has much smaller variance: 



SO that /f ~ 0.06 £2/5 ^ js 



s\A 6^ 



0.0039 £\ (108) 



2. The tensor part, 

The contribution of the tensor inhomogeneity to the statistical anisotropy, i^^, is a new 
result of the present paper. The relevant expression is given by formula (95) with To given 
by (82). One can write 

= {ns-l){ll + Jl) , (109) 
II = i^'n.Vg[T,-n(T,-T)]dT+{n^-n}, (110) 

Jl V{^h + -^h) [r, -n [n - r)] + {n ^ -n} . (Ill) 

The quantity is calculated by recalling definitions (63) with the result 

1 f^^ - 

Jl^-n^nn h'rir) in - t) dr . (112) 



It has only quadrupole dependence on the unit vector n. 
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The quantity is given by 

\ r dr I t^Mt) e.,(p,<T)nV(pn)sin[pn(Tfe-T)]e-^'/-V. + H.c. 

(113) 

It can further be expanded in polynomials in the unit vector n in a way similar to (105), 
(106): 

oo 
m=0 

where 



X / tSj^pIt) E e„(P,<')«V(pn)'<'"+"e-''/='"'ap,, + H.c. (116) 



The variance of the quadrupole term (112) can be estimated using (48), (53) and (67) as 
follows: 

(J^)) dnin - n) J dT,{r, - r,) dpp^ 

^ ^ / dnin - n) / dr.in - T,){a,H,YHl h'^. (ii6) 

This variance is independent of the effective cutoff parameter (cS in Eq- (67) due to the scale 
invariance of the spectrum of tensor perturbations. For our model (3), using the formulas 
of Sec. II A, we have 

^''^/-IbnJ^. a- low -lOW' 

where the last estimate takes into account the relation Hf — if ^ ~ Hf. 

In the same way, one can estimate the variances of the multipole terms (115): 

2\ 2e4(-+i)G .da 



mk,n)Y) ,/ H'- 

^ ^1 7r(4m + 5)(4m+ 7)(4m + 9)(m + l)2[(2m + l)!]^/a, a 

27r(4m + 5)(4m + 7)(4m + 9)(m+ l)2[(2m+ 1)!]^ ml ' ^ ' 
Again we see that the term with m = is dominating; for example. 



T\2 



0.0017 £\ (119) 
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so that ~ 0.04 £2/7^ < 

Quite similarly we can calculate the variances of the derivatives of the relevant quantities 
with respect to In/c, which turn out to be relatively very small: 

This estabhshes very weak dependence of v'^ on the absolute value of k for fixed n — 'k/k. 



C. Summary of the main results of this section 

In this section, we have found that the large-scale inhomogeneity (and anisotropy) of the 
inflationary universe due to the generation of quasiclassical scalar and tensor perturbations 
on super-Hubble spatial scales leads to statistically anisotropic corrections to the power 
spectrum of the post-inflationary primordial scalar perturbations. The correlation function 
of the relativistic potential $ has the form 



($(x)$(y)) = 




where i^k = > ^^<i the quantities and describe the influence of scalar and tensor 

quasiclassical inhomogeneities, respectively. The leading contributions to these quantities 
are given by 

= {ns - 1) (Agn^n^' + Kf^^^n'n^n'n') , (123) 
vl = (715 - 1) {^In'n^ + Kl^iu'n^n^v}) , (124) 

in which the quantities A;^;^ very weakly depend on the wave number k. One can introduce 
the combined quantities 

A,, = Ag + Ag , K^ki = A^jki + ^Jjki , (125) 

so that 

i^k = {ns - 1) {Aijn'n^ + Kijkin'nPn^v}) . (126) 

Results (123), (124) or (126) are proportional to {ns — 1), the deviation of the scalar 
spectral index from its value for the exactly scale-invariant form. 
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The quantities A^;-^ are given by the following expressions: 



Ag = 2 £ dr (r, - r) j A^^^^{r)p,p,e-^/^^'' + H.c. , (127) 

Ag,, = -\ £ dr in - r f j ^$^(r)p,p,.p,p,e-^/-^ap + H.c. , (128) 

AS = 2 / K{n-T)dT, 

^Iki = 2 / ^^(^'^-^) / 7^^f(^) E efe(p,<7)pfePoe~^^"V- + H-c., (130) 



(129) 

o-=+,x 

where, according to (67), the quantity h'lj is given by 

= ^ / E %(P' + H-c. (131) 

The variances of the anisotropic parts of the power spectrum are estimated as 

(AS"'"')') = {(Ilf) - (H^ - Hi) , (132) 

(A|„rfrf«'n')=) = {(irf) - 5g^G {Hf - Hi) , (133) 

(A>VW)>(W)>£^^1^. (135) 

The total anisotropy is dominated by the tensor contribution, which is clear from the 
ratio of their variances: 



(Agn^n^y) 2 ((Ag,,n^n%V)') 

^ ~ < 1 , -) -f ~ lO-'e^-l < 1 . (136) 



2 



Within the tensor (as well as within the scalar) contributions, the quadrupole component 
dominates: 

/Aj.H^VnV)'^ ^ ^4 (^(Ag;.;n^n%V)'^ ^ ^4 ^^^^^ 



(Agn^nO') ^21' ((AgnW)^) 259 

The full traces of the symmetric tensors Afj and Ajj^i are equal to zero, so that the quantity 
being averaged over the directions of n gives zero. On the contrary, the quantity gives 
a nonzero contribution to the isotropic part of the power spectrum: 

= ^ / ^k^^n = {ns - 1) Qa^' + ^Af;'^-) , (138) 

which, however, is rather small in view of the estimate (100). 
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VI. NUMERICAL ESTIMATES AND THE TRANS-PLANCKIAN ISSUE 

Observing the variances (132)-(135), one can see that they depend crucially on the value 
Hi of the Hubble parameter at the time r,, which is the moment of time where the initial 
conditions for quantum fluctuations are specified as S'k(rj,x) = Tk(Ti,x) = (see Sec. IV B). 
The evolution of the scalar modes is traced beginning from this moment of time roughly 
until the moment Tfc, which, we remember, is defined by the condition ak\Hk\^/'^ — k ; this is 
roughly the moment of the Hubble-radius crossing at the inflationary stage. 

The question is: at what time can one set the initial conditions? First of all, we re- 
member that our calculation is based on perturbation theory in terms of the quasiclassical 
inhomogeneities; this implies that we can use our expressions (132)-(135) without signifi- 
cant corrections as long as they are smaller than unity. The largest of them, (134), then 
determines the highest initial value for the Hubble parameter: 



^1 2 ^ 3vl0 2 /ior,\ 



where Lp^, is the typical value of the infiaton field that determines the boundary of the self- 
regenerating infiation, as described in Sec. IV A [see Eq. (58)]. As noted in that section, 
this is still much lower than the boundary of the Planckian energy density [sec (59)]. Thus, 
we might expect the variance (134) of the quadrupole part of the statistical anisotropy to 
be of any value between zero and unity, while that of the next leading term (135) to be an 
order of magnitude smaller. For ~ variances (132) and (133) of the contributions 
from scalar inhomogeneities are small and can be neglected:^ 

The dispersion of the anisotropic part in this case is due to the contribution from tensor 
inhomogeneity and is estimated as a^^^ ~ {ns — 1) ~ 0.04. We emphasize, however, that 
this is only the value at which our perturbative approach will possibly require corrections; 
in reality, if 0j ^ 0*, the anisotropy may be larger. 

There is another important issue that affects the resulting statistical anisotropy. The 
point is that the physical wavelength X — a/k oi a particular mode is much smaller than 
the Planck length ip ~ G^/^ during most part of the inflationary universe. This situation is 
illustrated in Fig. 1. A natural question arises whether it is legitimate to use fleld theory in 



^ The scalar and tensor contributions to the anisotropy become equal by the order of magnitude if the 
initial value of the Hubble parameter is as low as Hi ~ mip, which is very close to the end of inflation. 
Thus, in all interesting cases, the dominating contribution to the statistical anisotropy comes from tensor 
inhomogeneities . 
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FIG. 1: Illustration of the trans-Planckian issue. The physical wavelength A = a/k of a particular 
mode is below the Planck value l-p — G^/^ during most part of the inflationary universe. The 
Planck-radius crossing \ = l-p occurs only about 14 e-foldings prior to the Hubble-radius crossing 
A = during inflation. 

cosmology on spatial scales below the Planck length. This is the so-called trans-Planckian 
issue. In the context of the origin of statistical anisotropy, it was pointed out in our papers 
[4], and recently it was a subject of discussion in several works in connection with the possible 
modification of the (statistically isotropic) primordial power spectrum [8] and cosmological 
particle creation [9]. 

To see how this consideration affects our results, let us calculate the values of the infia- 
ton field ipk (at the moment of Hubble-radius crossing, where, according to our definition, 
k = a\H\^^'^) and (fp^ (at the moment where the physical wave number crosses the Planck 
boundary, k ^ aG^^^"^) for our quadratic inflaton potential. Using the formulas of Sec. II A, 
we have 

47rG'(^2 ^21n^^^ ^ 100, ATrGipl,, = AnGifl + In ^ 128 . (141) 
k Gmfp 

Together with (34), these estimates imply that the Planck-radius crossing occurs only about 
2ts:G {^Ppk — ^p1) ~ 14 e-foldings prior to the Hubble-radius crossing during inflation. The 
difference between cppk and is rather small, and, if we set the isotropic initial conditions 
for the modes at the Planck-radius crossing, so that ipi = ippk, then, for variance (134) of 
the leading contribution to the anisotropy, we will have 

\^ ■' 'I lOvrm:! ^ 
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which is extremely small. The statistical anisotropy in this case is negligible: (f^) ~ 10~^^. 

There are several arguments in favor of the possibility of working in the trans-Planckian 
frequency region. First of all, one can note that the wavelength and frequency of a particular 
mode are not Lorentz-invariant quantities and can locally be changed to any value by the 
transformation of coordinates. Thus, the authors of [9] argue that nothing dangerous results 
from this procedure in a locally Lorentz-invariant theory as long as the invariant k'^k^ <^ Mp. 
It is also pointed out there that a similar trans-Planckian problem arises formally in the 
calculation of electron-positron pair creation in a constant electric field Eq if the gauge 
A = — Eoi is used [23]; in this case, in order to obtain correct results for the effective 
action, it is even obligatory to consider field modes with arbitrarily high momenta. Similar 
arguments in cosmological context are given in [24], and the authors of [25] consider the 
trans-Planckian effects in quantum gravity and argue in favor of the asymptotic safety in 
this region. 

The trans-Planckian issue can be put in a slightly different way. Given that there exists 
a formal moment of Planck-radius crossing for every mode (see Fig. 1), one can ask whether 
the quantum state in which a particular mode finds itself at this moment of time will be 
sensitive to the past history of the inflationary expansion. The unknown physics taking place 
in the trans-Planckian region might have lead to an effect similar to that of our field-theory 
consideration, producing the modes Xk and at the Planck-radius crossing similar to those 
calculated in this paper. This issue will hopefully be clarified in a reliable quantum theory 
capable of dealing with sub-Planckian spatial scales. 

One can note that quite a similar trans-Planckian problem arises in the theory of quantum 
radiation from a black hole which is formed as a result of gravitational collapse: since 
all outgoing modes originate as incoming modes that arrive at the black hole before the 
formation of the event horizon, the frequency of the incoming modes diverges as the time 
of the corresponding outgoing modes goes to infinity (see [26] for a review). While the 
prediction of quantum thermal radiation of black holes hardly can be tested experimentally 
because of its extremely small value for black holes with realistic masses, we have a different 
situation in the case of the cosmological prediction made in this paper, which, in principle, 
can be verified by observations. We discuss such testable predictions in the next section. 

VII. OBSERVABLE EFFECTS 

Statistical anisotropy of the primordial power spectrum will broadly be manifest in the 
large-scale structure and cosmic microwave background. The first effect was discussed in our 
previous paper [4], while the statistical anisotropy of the CMB was a subject of a number 
of recent studies [10-16]. 
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A. Statistical anisotropy of the LSS 

Since the leading contribution to the statistical anisotropy is of quadrupole form [the 
corresponding variance is given by (134)], the formulas of [4] are directly applicable to the 
present case. There is a linear homogeneous and isotropic connection between perturbations 
in the relativistic potential $ and density perturbations 6 = Sp/p on sub-Hubble spatial 
scales. Therefore, the anisotropic factor z/k is inherited in the power spectrum of the density 
contrast. The corresponding correlation function is given by 

e(x) = eo(a;)+ei(x), (143) 

where 

«.)=/|^^^e-^^/d..^.^!5^ (144) 

and 

6(x) = / (03^^-^''^ (145) 

are the isotropic and anisotropic parts, respectively. The anisotropic quantity i/k is domi- 
nated by the first term in the brackets of (126) in which Aj^ Kj-. In view of the fact that 
the product {ns — 1) Aj^ depends very weakly on the wave number one can treat it as a 
constant and take out of the integral in (145). Since the matrix A^ is traceless, one gets the 
following expression for the anisotropic part of the correlation function: 

i^{-s.)^{ns-l)lyiif{x)x"x\ (146) 

where x = x/x and 

3 f n kx COS kx sin kx ,„\ 

fix) = eo(^) + dkkX . (147) 

The anisotropic part of the correlation function has the specific form (146), (147) which, 
in principle, can be tested by observations. Measurement of the quadrupole of the power 
spectrum can simply be done by comparing the amplitudes of Fourier modes in different 
directions. The standard error to the quadrupole coefficients qij — {ns — 1) using the 
SDSS data is estimated in [12] to be aq ~ 10"^. 

The velocities on large spatial scales are connected with the relativistic potential as 
V oc V$. Therefore, one can also test the statistical anisotropy by using the velocity 
correlation function (Tij(x) = (xq — |) vj (xq + f ))■ Of a particularly simple form is the 
dispersion of the velocity components at one point: 

a,, (0) = {Vi (xo) Vj (xo)) = ^ (v^) (s,j + ^K,^ . (148) 
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This relation is valid on large scales, where the evolution of density perturbations is still in 
the linear regime, and, in principle, can be used to determine the matrix Aj^. 

The effect of nonlinear evolution of initially anisotropic spectrum was recently under 
investigation in [27] ; the authors reported the suppression of primordial anisotropy by ~ 7% 
in quasilinear regime on small spatial scales, k ~ 0.1 Mpc~^. 



Statistical anisotropy of the primordial power spectrum in the form (126) falls into the 
category that was recently studied in detail in [12]. It is given by the general formula 



where is the isotropic part of the anisotropic power spectrum Pk of the density profile 
S — 5p/p, and gLuik) are the coefficients of the expansion of the statistical anisotropy into 
spherical harmonics Ylm{^), where, as usual, n = k//c. In our case, the sum in (149) is 
dominated by harmonics with L — 2, with coefficients g2M ~ {ns — 1) \j that depend on k 
very weakly. 

The primordial inhomogeneity described by the power spectrum (149) leads to the 
anisotropic distribution of the CMB temperature T{n) on the sky; the relevant relation 
can be expressed as follows: 



where n denotes a direction on the sky, and Qi{k) is the kernel of the linear relation between 
T(n) and Sk describing the contribution to the Ith temperature moment from the wave 
vector k [with these conventions, Qi{k) is real]. 

After expansion of the temperature map T'(n) into spherical harmonics 



B. Statistical anisotropy of the CMB 




(149) 




(150) 




(151) 



Im 



one can calculate the expectation value (covariance matrix) 





LM 



with 




(153) 
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/ dkePk9LM{k)Qi{k)Q,{k), 



(154) 



D 




(155) 



The quantities D^, are the generahzation of the coefficients Q to the case of statistical 
anisotropy. One can note that, in the case where the dependence of QLuik) on k can be 
neglected, there arises the relation Df(F — QlmWu' with Ww depending only on the isotropic 
cosmology (in particular, Wu — Ci). Thus, the combinations D^M /Ww can be regarded as 
convenient estimators for the coefficients qlm in this case. The estimators for Df(/^ were 
constructed in [12] and essentially coincide with the bipolar-spherical-harmonic coefficients 
of [10]. For the case L = 2 of interest for our work, it is shown in [12] that the variances 
with which the quantities g2M can be measured using the temperature maps of WMAP 
and Planck are ag ~ 1.2 x 10~^ and 3.8 x 10~^, respectively. This is comparable to the 
highest theoretical estimates that we get within the limits of our approximation, linear in the 
background inhomogeneity of the inflationary universe. In a complete nonlinear calculation, 
therefore, the resulting statistical anisotropy may well exceed the current threshold of its 
detection. 

Our theory does not predict any significant dipole statistical asymmetry of the CMB or 
LSS, which is a simple consequence of the fact that the correlator ($(x)$(y)) in our approx- 
imation is a symmetric function of the difference x — y. In fact, statistical inhomogeneity in 
our theory would be obtained in the next order of expansion of the phases in (92) in powers 
of X and y, and one can show that this effect is very small. Specifically, the spatial variance 
of the quantity can be estimated as ((x • Vukf) ~ {ns - if {kxfGml ~ lQ-^^{kxf. 
Even if one takes the comoving length x as large as the size of the observable universe, 
X ~ 3000 Mpc, and comoving scale k as large as /c ~ IMpc"^, one gets the estimate 
((x • Vi^k)^) ~ 10~^. In particular, this imphes that our theory by itself does not predict 
any significant dipole modulation of the CMB power spectrum, which is suggested to ex- 
plain the observational evidence of the hemispherical power asymmetry in the WMAP data 
[13] (see [14] for a recent discussion of this issue). However, the quadrupole modulation 
of the CMB sky, which is present in our theory, together with anisotropy in CMB data 
(caused, e.g., by the cut-sky masks used to exclude foregrounds, see [15]), might be respon- 
sible for this effect. Indeed, in this case, effectively nonantipodal regions of the sky may be 
involved in the comparison of the power spectra of two antipodal hemispheres, and then the 
quadrupole asymmetry may produce the effect of antipodal asymmetry. In our case, the 
effect is described by an almost scale- independent traceless matrix Aj^ , hence, it is expected 
to exhibit universal axes of asymmetry for statistical anisotropy on all angular scales. This 
question remains to be investigated in all detail. 
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VIII. SUMMARY 



In this paper, we considered the evolution of small-scale quantum fluctuations on the 
background of an inhomogeneous inflationary universe and calculated the resulting primor- 
dial power spectrum. The background metric inhomogeneities during inflation are generated 
in a self-consistent way on super-Hubble spatial scales. 

We have seen that the resulting power spectrum of scalar primordial perturbations ac- 
quires the statistically anisotropic part and can be put in the form (1), (2). The quantities 
Ajj and Aj^fc; in (2) very weakly depend on the wave number k. Being a result of metric 
inhomogeneities generated during inflation, they represent Gaussian random variables with 
zero average; in our universe, we observe only one particular realization of these quantities, 
hence, we can only speak about their statistical properties. 

The variance (which characterizes the expected magnitude) of the statistical anisotropy 
crucially depends on the moment of time at which one specifies the initial conditions for 
the modes describing the vacuum state. It is here that we encounter the trans-Planckian 
issue, which we discussed in Sec. VI. If one specifies locally homogeneous and isotropic 
initial conditions for quantum modes at the moment of Planck-radius crossing, then the 
statistical anisotropy which we are talking about turns out to be very small, with variance 
{u^) ~ 10"^"^. However, one can argue in favor of the possibility of working in the trans- 
Planckian domain of wave numbers (assuming exact local Lorentz invariance) and specify 
the initial conditions much earlier along the inflationary epoch. In this case, the expected 
value of the statistical anisotropy can be appreciable; our perturbative approach breaks 
down at the values (u^) ~ 10"^; in reality, the anisotropy might be considerably higher and 
might be detectable. 

Where the statistical anisotropy is sufficiently high, it is mainly caused by the influence 
from the background metric inhomogeneities of tensor type. In this case, the full traces of 
the tensors Ay and Aij^i in (2) are equal to zero, and the ratio of the corresponding variances 
is estimated in (137) to be 



The speciflc scale- invariant form (2), dominated by the quadrupole term, might facilitate 
the search for the predicted statistical anisotropy in the large-scale structure and cosmic 
microwave background. Its detection may be regarded as another nontrivial conflrmation of 
the chaotic inflation scenario and elucidate some aspects of the trans-Planckian issue. 




(156) 
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APPENDIX A: DERIVATION OF EQ. (93) 

Since the contributions from S'k and Tk add up in the linear approximation, we illustrate 
the derivation of (93) from (92) by considering only one of these quantities. Consider then 
the integral 

I = J ^^^|^j.|2e'^(''-y)-'^k('"fc'^)+^^k(^'^'y) . (Al) 

According to the procedure described in Sec. VA, we expand the exponent into a series 
in the spatial coordinates retaining the leading terms: 

fk (Tk, x) ^ kfo (Tk, 0) + A;x • Vfo {Tk, 0) + ifi (r^, 0) . (A2) 

The exponent of (Al) becomes 

ik(x - y) - ifk (Tfe, x) + ij;: (Tfe, y) ^ i(k - (5k) (x - y) + 2f, (r^, 0) 

= i(k - 5k)(x - y) + 2T, (Tfe, 0) + 2n • VTo (r^, 0) + 2Q (r^k, 0) , (A3) 

where 

<5k = A;VTo(rfc,0) , (A4) 

and we have used (88). In what follows, we omit the argument (r^, 0) which is going to be 
present in all functions. The integral (Al) then becomes, in the linear approximation in the 
quasiclassical field h^j, 

J^IM'' (1 + 2Ti + 2n • VTo + 2Q) e'^-'^-^'^^^^-y^ . (A5) 

Now we change the integration variables from k to k — (5k. We have to take into ac- 
count the transformation of the amplitude \Ak\ and the Jacobian of transformation to new 
variables, all in linear approximation with respect to 5k defined in (A4). The first part gives 



\Ak+5k\' = \Ak\' + -^Sk = \Ak\' ( 1 + • VTo ) . (A6) 



The Jacobian in linear approximation is calculated using (79), (82) and (83) as follows: 
J = l+tr-^ = l + n.VTo-2Ti-g + [n • VQ - (n • V)' To] ^ . ( 
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Combining (A5), (A6) and (A7), we obtain 



I = 



where 



d\n{k^Ak\ 



- " kT' - n . Vro + Q + [n . VQ - (n . V)' r.J . (A9) 

Considering the second integral in (92), we arrive at an expression similar to (A8) but 
with //^j. replacing /i^. By adding these two integrals with the factors 1/2, we obtain 

'A,mi + ul)^^(-y\ (AlO) 



dTk 



(2vr) 



where 



I (/^k + /^-k) = g + J (ns - 1) n • V [(To)k - (To)_k 



2 

^n.V)n(T„V+(T.)_,l-i|^. (All) 



2 

The variance of the local quantity Q can be estimated using (66) and (53), and we get 



<^'>" / ^s\hAn)f^GHlln'-c.G^my,c.lO-'\ (A12) 



Because of its smallness, the local term Q can be neglected in (All). 

The last term in (All), containing the derivative drk/dlnk ~ \Hk\^^'^/kHk, arises from 
the Jacobian (A7) as a result of the dependence of To on the time Tfe. This term is also much 
smaller than the second, leading term in (All), which should be clear from the fact that it 
is of higher order in the expansion parameter {k£)~^. One can calculate the variance of this 
term to find that it is of the order Gm^ ~ 10~^^. Thus, we should drop this term as well 
from the correction to the power spectrum. Eventually, we get our result (95). 

Quite similarly, one obtains contribution (94) from the phase S^. The local quantity 
which is analogous to Q in (All), as well as the expression similar to the last term in (All) 
are both much smaller than the leading term (94). 

In principle, expression (93) for the correlation function should be revised in view of the 
fact that the spatial coordinates x entering the metric at the end of inflation are not exactly 
Euclidean because of the presence of the large-scale quasiclassical inhomogencitics in metric 
(64). However, proceeding to the Euclidean coordinates will generate additional local terms 
of the form and order of $c and Q in the anisotropic part of the correlation function, which 
thus also can be neglected. 
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APPENDIX B: ESTIMATES OF THE VARIANCES 



The variance of integral (97) is given by 



(7^)') = 4 Tdn rdT2 

X J ^^$p(ri)$;(r2)(pn)2 sin [pn {n - n)] sin [pn {n - Ta)] e'^i^^^) , (Bl) 

where ai = a(Ti), Hi — H{ti) etc. The cutoff value of momentum integration is sufficiently 
low, so that the argument of the sine in the effective integration region satisfies 

|pn(rfe - r) I < eaH{Tk - r) < 1 . (B2) 

Taking into account (29) and (31), we can estimate variance (Bl) as follows: 

dn^^ in - n) £ dr,^, (r, - r^) / ^ \A/ (pn)VKi^-^^) 



57r2 



5 



I" dn^ (Tfe - n) r dT2 in - T2) 4hI , (B3) 



where, in the last equality, we have used equations (6) and (7). Now using the approximate 
relation ^ ^ 

^ aH afei^fc ' ^ ^ 

and taking into account the slow variation of the Hubble parameter as a function of the 
scale factor a during inflation, we get the estimate 

,,cx2\ 4£4G2 p Ae^C^ .2rr, "^^^G^ f"^' dV iif) , 

<^^ ' ) " — 1 ^ v — Z ^ - -IT i 

= ^K-V,)c.^G{H!-Hi). (B6) 
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